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COMPUTATIONAL ENGINE STRUCTURAL ANALYSIS 


Christos C. Chamls and Robert H. Johns 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland , Ohio 44135, U.S.A. 


SUMMARY 

A significant research activity at the NASA Lewis Research Center Is the 
computational simulation of complex multidisciplinary engine structural prob- 
lems. This simulation Is performed using computational engine structural 
analysis (CESA) which consists of Integrated multidisciplinary computer codes 
In conjunction with computer post-processing for "problem-specific" applica- 
tion. A variety of the computational simulations of specific cases are 
described In some detail In this paper. These case studies Include (1) aero- 
elastic behavior of bladed rotors, (2) high velocity Impact of fan blades, (3) 
blade-loss transient response, (4) rotor/stator/squeeze-f llm/bearlng Interac- 
tion, (5) blade-fragment/rotor-burst containment, and (6) structural behavior 
of advanced swept turboprops. These representative case studies were selected 
to demonstrate the breadth of the problems analyzed and the role of the com- 
puter Including post-processing and graphical display of voluminous output 
data. 


INTRODUCTION 

The need for lightweight, durable, fuel-efficient, cost-competitive air- 
craft requires engine structures which are made from advanced materials 
Including composites, and which resist higher temperatures, maintain tighter 
clearances and reduce maintenance costs. These requirements tend to Increase 
the complexity of aircraft engines which are already Intricate structural 
systems consisting of dynamically Interacting components. Prediction of the 
Integrated response of these systems requires complex analytical models (finite 
element, for example) as well as complex static, periodic, transient, and 
Impact analyses. A major engine structures program (fig. 1) Is currently being 
conducted by NASA Lewis, with emphasis on computational engine structural 
analysis (CESA). Several Independent computer codes for specific aspects of 
engine structures modeling and/or structural analysis are being developed as 
part of this program. These codes Include: (1) Interactive finite elements 

for bearing/shaft Interaction, (2) engine structural component modeling, (3) 
automatic finite element generators, (4) nonlinear analysis with various levels 
of sophistication, (5) nonlinear constitutive relationships, (6) automatic 
thermal load transfer, (7) blade-loss transient analysis, (8) aeroelastlc 
analysis of bladed disks, (9) multimode Impact analysis, and (10) structural 
tailoring of engine blades. 

Several of these computer codes were briefly described In a previous paper 
(ref. 1). The objective of this paper Is to Illustrate typical computational 
results obtained from these codes for: (1) aeroelastlc behavior, (2) high- 

velocity Impact, (3) blade-loss transients, (4) rotor/stator/squeeze fllm/- 
bearlng Interaction, (5) blade fragment containment, and (6) advanced turbo- 
prop structural behavior. The role of the computer In the solution of these 



complex multidisciplinary problems, the post-processing, and the graphic 
display of voluminous results obtained therefrom are emphasized. 


AEROELASTICITY 

The multidisciplinary analysis required for predicting aeroelastic 
behavior including flutter is performed by a special-purpose code (ref. 2). 
This code couples appropriate aerodynamics with the structural response of 
bladed, shrouded disks. Flutter is an Important aeroelastic design consider- 
ation for fan blades. When it occurs it subjects a blade to high-amplitude 
cyclic loading which can rapidly degrade its structural integrity and lead to 
blade failure. 

Blade failure Induces rotor imbalance which, if of sufficient magnitude 
can cause rotor stage failure and/or bearing failures, both of which can lead 
to catastrophic engine failure, if only very Infrequently. 

The computer code predicts the flutter speed and the blade motion near or 
at this speed. The code uses cyclic symmetry substructuring concepts and the 
widely-used NASTRAN structural analysis computer code. A computer plot of a 
typical sector of a bladed shrouded disk analyzed is shown in figure 2. Com- 
puted graphical results of the motion of this sector near flutter speed are 
shown exaggerated in figure 3. The character of the motion and location of 
high bending (fatigue) stresses are readily observable in this figure. 


HIGH VELOCITY IMPACT OF FAN BLADES 

High-velocity Impact which can result in foreign object damage (FOD) is an 
Important design consideration in aircraft turbine engine fan blades. Some 
research activities at NASA Lewis focused on developing technology for impact- 
resistant fan blades using advanced design concepts and several computer codes 
(ref. 1). 


A computer-produced three-dimensional finite element model for one of these 
advanced blade concepts is shown in figure 4. Stress contour plots obtained 
using this finite element model are shown in figure 5 for several blade depth 
locations or "surfaces." Impact location and blade root stresses obtained 
using modal synthesis for transient analysis are summarized in figure 6. 

The amount of information associated with high-velocity Impact analysis is 
very voluminous. Computer post-processing and computer graphics including 
computer-produced movies is the practical way to handle it. Furthermore, 
computer color graphics enhance the interpretation of the output relative to 
displacement magnitudes, stress distributions and magnitudes, and transient 
motion of the Impacted blade. 


BLADE-LOSS TRANSIENTS 

Partial or full blade loss due to Impact, for example, causes stage 
tmbalance which can Induce rotor transient whirling. This type of imbalance 
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and the resulting transient and steady-state motion must be predictable so 
that appropriate tolerances and operating restrictions can be provided. Pre- 
diction of blade loss transient response requires an integrated structural 
dynamic analysis where the bladed rotor, shaft, and bearings are all modeled 
in the analysis. A computer code (Turbine Engine Transient Response Analysis 
[TETRA] ref. 1) was developed for such an Integrated analysis. 

TETRA uses the component element analysis method with progressively reduced 
model complexity as shown in figure 7. TETRA output is also voluminous. Com- 
puter graphics is the practical way to display this type of output. Typically 
TETRA produces displacement/time history plots as shown in figure 8, or cor- 
responding frequency-domain plots as shown in figure 9. Both of these plots 
are valuable and necessary in interpreting the blade loss transient response 
and in assessing the adequacy of the design. 


ROTOR/STATOR/SQUEEZE-FILM/BEARING INTERACTION 

Squeeze-film damper bearings are designed for specified blade-loss transi- 
ents and rotor imbalance and vibration. Rotor/squeeze-film/bearing/ stator 
interactions are inherently nonlinear because of: (1) the inherent behavior 

characteristics of squeeze films, particularly during transient situations, 
and (2) potential structural interactions wherein either large deformation 
kinematics or material nonlinearity (plasticity) are excited. 

Because of such behavior, the overall rotor-bearing-stator computational 
simulation must be able to incorporate the various sources of nonlinearity. 
Additionally, to enable efficient solutions in situations wherein the behavior 
of the model components is linear, the overall simulation scheme must incorpo- 
rate substructuring capabilities. Furthermore, since transient problems need 
to be considered, such features must be accommodated by the various integration 
algorithms used to solve the governing model equations. 

A general purpose squeeze-film-damper/interactive-force finite element 
model has been developed as part of the NASA Lewis program (refs. 3 to 4). 

This finite element model was implanted in a general purpose finite element 
computer code. Three different numerical integration methods were also incor- 
porated to solve the governing structural dynamics equations of the interactive 
finite element. The general purpose code is used subsequently to predict the 
structural dynamic response of the rotor/stator coupled structure subjected to 
imbalance and impulse type excitations. The structural dynamic responses pre- 
dicted include: (1) bearing/rotor trajectories, (2) stator trajectory, (3) 

rotor orbit, and (4) force, velocity and acceleration histories at a given 
location in the coupled structure. Three-dimensional post-processors are used 
to graphically display the predicted results in isometric views. Typical com- 
puter post-processing results are shown in figure 10 for a shaft/single bearing 
interaction, and in figure 11 for a shaft/multibearing interaction. These 
graphical results clearly indicate squeeze film pressures, orbit motion, orbit 
stability, clearances, and respective margins. 
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BLADE-FRAGMENT/ROTOR-BURST CONTAINMENT 


Blade-fragment containment and rotor-burst containment are critical design 
requirements for turbine engines for commercial aircraft. The computational 
simulation of blade-fragment and/or rotor-burst containment require complex 
analyses Including fragment tracking, nonlinear material behavior, and predic- 
tion of finite strains, large rotations and large displacements. These complex 
analyses were developed as part of the NASA Lewis program and Incorporated 
Into an Integrated computer code using nonlinear finite element theory (ref. 5) 

The containment computational simulation using this code starts with the 
Initial configuration of the containment ring, rotor speed, fragment size, and 
fragment location. The simulation proceeds to determine the fragment kinetic 
energy, momentum, direction, and Impact point on the containment ring. The 
transient structural response of the ring from this contact point Is tracked 
as a function of time. The transient structural response generally Includes 
displacements, strains and stresses. 

Because of the Importance of fragment containment, extensive experimental 
studies have been conducted on specific rotors and containment rings. The 
computational simulation results for one of these cases are Included herein. 

The containment ring geometry and properties and the rotor characteristics are 
shown In figure 12. A computer plot of the shape of the deformed ring at about 
Its maximum deformation Is shown In figure 13. Computer plots for the stresses 
at the corresponding time are shown In figure 14 for the outer surface and In 
figure 15 for the Inner surface. 


ADVANCED PROPELLER STRUCTURAL BEHAVIOR 

The use of thin, highly swept, twisted propeller blades Is Important to 
the development of fuel-efficient and quiet advanced aircraft. High sweep 
angles produce significant reductions In noise and are therefore a desirable 
design feature. However, blades of this type (thin, highly swept, twisted) 
exhibit complex structural response under a centrifugal force field, which 
requires special analysis techniques for accurate characterization. These 
techniques are required to establish the structural dynamic response of the 
turboprop blades at operating rotor speeds Including the avoidance of possible 
leading edge buckling under centrifugal loads. Computational simulation 
studies (refs. 6 to 8) were performed at NASA Lewis to determine analytically 
the structural behavior and possible leading edge buckling of advanced, highly 
swept turboprop airfoils In centrifugal force fields. Computational simulation 
studies were also performed In order to Identify any advantages of using "high 
performance" composites In order to change the regions of Instability. The 
theoretical studies were performed using an In-house program (COBSTRAN, ref. 1) 
designed for composite blade analysis. Schematics of the rotor, the propeller 
and the finite element model used In the analyses are shown In figure 16. 

Computer output results for rotor speed effects on the vibration mode 
shape and potential for leading edge buckling are shown In figure 17. Computer 
plotted mode shapes of the propeller blade superimposed on a Campbell Diagram 
are shown In figure 18. Display of the amount of computer output Information 
summarized In figure 18 would be Impractical without computer post-processing 
and graphical display. The availability of all this Information In one figure 
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makes It convenient to assess the adequacy of the design on a quantitative 
basis . 


SUMMARY OF RESULTS 

The complex behavior of Integrated aircraft engine structures requires the 
solution of interacting multidisciplinary problems. Some recent research 
activities at NASA Lewis focus on computational engine structural analysis 
(CESA). CESA Is essentially the computational simulation of complex engine 
structural behavior using multidisciplinary Integrated computer programs In 
conjunction with extensive computer post-processing and graphical display of 
computer output. Representative case studies were selected and are described 
In some detail to Illustrate the role of the computer on the computational 
simulation of complex multidisciplinary engine structural problems. These 
case studies Include: (1) aeroelastic behavior of bladed rotors, (2) high 

velocity Impact of fan blades, (3) blade-loss transient response, (4) rotor/ 
stator/squeeze-f llm/bearing interaction, (5) blade-fragment/rotor-burst con- 
tainment, and (6) structural behavior of advanced swept turboprops. These 
representative case studies further demonstrate the Indispensable utility of 
computer post-processing and graphical display of voluminous analysis results 
which make It possible to readily assess the adequacy of a design on a quanti- 
tative basis. 
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Figure 1. - Engine structural components amenable to computational simulation. 
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Figure 3. - Aeroelastic finite element analysis of bladed shrouded disks 
(bell code). 
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Figure 4. - Superhybrid composite blade finite element model. 












Figure 9. - Demonstrator model, 100 g-in. sudden unbalance at 
the fan at 12 000 rpm (spectrum-log mag. versus linear fre- 
quency). 
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Figure 10. - Bearing damper interactive element - pressure and orbit responses. 
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Figure 11. - Structural dynamics analysis using rotor/stator/ squeeze-film/bearing interactive finite element. 
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Figure 12. - Geometric, test, and modeling data for the 4130 steel containment ring 
subjected to tri-hub T58 rotor burst in NAPTC test 201. 
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Figure 13. - Deformed configuration 1180 |isec after burst, 
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surface strain 1180 Msec after burst. 




TIP 





744 ELEMENTS 

Figure 16. - Turboprop stage propeller and finite element model. 
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Figure 17. - Structural behavior of advanced turboprops. 
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Composite turboprop vibration analysis. 
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